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Supplementing the Liquid Drop Model (LDM) Hamiltonian, written in the intrinsic reference 
frame, with a sextic oscillator plus a centrifugal term in the variable /3 and a potential in 7 with a 
minimum in ~, the Schodinger equation is separated for the two variables which results in having 
a new description for the triaxial nuclei, called Sextic and Mathieu Approach (SMA). SMA is 
applied for two non- axial nuclei, 180 Hf and 182 W and results are compared with those yielded by the 
Coherent State Model (CSM). As the main result of this paper we derive analytically the equations 
characterizing SMA from a semi-classical treatment of the CSM Hamiltonian. In this manner the 
potentials in /3 and 7 variables respectively, show up in a quite natural way which contrasts their 
ad-hoc choice when SMA emerges from LDM. 



I. INTRODUCTION 

Many properties of the low lying spectra of even-even 
nuclei can be described in terms of specific dynamic sym- 
metries [l|, associated with a definite nuclear phase. 
The transition from one symmetry to another is therefore 
interpreted as a phase transition. On the path of such a 
transition a critical point is met. The spectroscopic prop- 
erties of the nuclei corresponding to the critical points are 
in general difficult to be described. 

In Ref. [3, 4j, it has been proved that on the transi- 
tion from the U(5) to the 0(6) symmetry a critical point 
exists for a second order phase transition while the tran- 
sition from the U(5) to the SU(3) symmetry has the fea- 
tures of a first order phase transition. In Ref.[5| it was 
proved that most of nuclei are mapped not on the border 
of the symmetry triangle introduced by Casten § but in 
the interior of the triangle. Examples of such nuclei are 
the Os isotopes Q. 

Recently, Iachello [8, 9] pointed out that the mentioned 
critical points correspond to distinct symmetries, namely 
E(5) and A (5), respectively. For the critical value of an 
ordering parameter, energies are given by the zeros of a 
Bessel function of half integer and irrational indices, re- 
spectively [Tol - [T2j . In Ref.[13| the X(5) description was 
extended to the first octupole vibrational band in nuclei 
close to axial symmetry and also close to the critical point 
of the U(5) to SU(3) phase transition. Another symme- 
tries, called Y(5) and Z(5), have been pointed out in 
Refs.fbH. [l5l|. The former symmetry corresponds to the 
critical point of the transition from axial to triaxial nu- 
clei while the latter one is related to the critical point of 
the transition from prolate to oblate through a triaxial 
shape. 

The nice feature of the critical point symmetry is that 
the description in the intrinsic frame is performed by two 
separated differential equations for beta and gamma de- 
grees of freedom. These equations are solvable and the 
solutions are irreducible representations for the specific 
symmetry. Moreover, apart from an overall scaling pa- 



rameter the energies are parameter free quantities. Since 
the idea of symmetries associated to the critical points 
of various phase transitions showed up, many attempts 
have been made to describe the two dynamic deforma- 
tions by solvable and separable differential equations with 
specific beta and gamma potentials. Since the triaxial 
nuclei might be considered as critical points for a phase 
transition from prolate to oblate shape one expects that 
they can be described by specific solvable models. Thus, 
a description of soft gamma nuclei around 7 — f with 
an oscillator potential in 7 and a Kratzer potential in j3 
has been developed in Rcf s . fi~6l4T8l j . 

Another model was proposed by two of us (A. A.R. and 
P.B.) in a previous publication [L9||. Indeed, the Liquid 
Drop Model (LDM) Hamiltonian [20| written in the in- 
trinsic frame is separated into two terms describing the /3 
and 7 variables. The potential in ft consists in a centrifu- 
gal term and a sextic oscillator potential, while the dif- 
ferential equation for 7 is that for the Mathieu function. 
Due to this feature we called the formalism developed 
there as the Sextic and Mathieu Approach (SMA). The 
proposed model was applied for five nuclei which exhibit 
the signature for triaxiality, 188 0s, 190 Os, 192 0s, 228 Th, 

Here we continue the study of SMA by describing an- 
other two nonaxial isotopes, 180 Hf and 182 W. In Ref. fl9| 
we remarked that the agreements with experimental data 
provided by SMA and the Coherent State Model (CSM) 
respectively, are of similar quality. Here we attempt to 
answer the question is there any reason for that to hap- 
pen? 

The objectives of this paper are described according 
to the following plan. In Section II, a brief review of the 
main ingredients defining SMA is presented. Section III 
is devoted to the CSM approach. It is shown, in Section 
IV, that SMA's equations are obtainable from CSM un- 
der certain circumstances. The numerical application is 
discussed in Section V, while the conclusions are drawn 
in Section VI. 



2 



II. SEXTIC AND MATHIEU APPROACH 

Adding to the Bohr-Mottelson Hamiltonian written in 
the intrinsic frame of reference [2(| a potential of a con- 
venient form: 



V(p,j) = V 1 (p) + ^V 2 (j) 

the equations for the variables P and 7 are separated 
from each other. Thus, the equation for /3 reads: 



P i dp H dp 



L(L + 1) 



f(P)=e f} f(p), 



The remaining terms, depend on 7 but also on P by 
means of the factor 1//3 2 . In order that the variable sep- 
aration is achieved, the mentioned factor is replaced by 
an average value l/(/? 2 ). Actually, in our concrete cal- 
culation this is considered to be a free parameter. The 
resulting equation in gamma variable, is: 
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where the following notations are used: 

= V2 ^ = 
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In order to solve the separated equations in f3 and 7 re- 
spectively, we have to specify the potentials V\{p) and 
^2(7). Vi(f3) is a sextic oscillator potential in /? and the 
corresponding differential equation is quasi-exactly solv- 
able. Indeed, by changing the function f(p) = P~ 2 tp{P) 
in Eq. (|2.2p one obtains: 
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One can check that this equation can be identified to 
the Schrodinger equation 

H x ip(x) = Eip{x) (2.6) 

associated to a sextic oscillator with a centrifugal barrier 
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The eigenfunctions provided by Eq. (|2.5p are written in 
the form: 
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(2.1) where N nj3t L are normalization constants and L (P 2 ) 

are polynomials of degree rip in p 2 , whose coefficients 
form an (M+l)-vector which satisfies an eigenvalue equa- 
tion, the corresponding eigenvalue being denoted by 



X^\L). Using the notations from Eq. (|2.4p . one obtains 
for the eigenvalues the following expression: 



E p (n Pt L) = — [Abs(L) + X^p(L) + < 
n =0,1,2, ...,M. 



(2.10) 



Here Uq denotes two constants to be fixed such that for 
the minima {Pmin > 0) °f * ne potentials vf (/3) and (/?) 
given by 

vl(p) = (b 2 - 4ac n )p 2 + 2abp i + a 2 p 6 + (tt = ±), 

(2.11) 

to have the same energy. Details about how to solve the 
eigenvalue equation for a sextic oscillator plus a centrifu- 
gal term can be found in Ref . [2TI] . 

Concerning the equation for the variable 7, Eq. (l2.3p 
can be reduced to the Mathieu equation [22J]. First we 
change the function 
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The equation for the new function is: 
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(2.13) 
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where L denotes the angular momentum and R is its 
projection on the axes OX. The potential in 7 is chosen 
to exhibit a minimum at 70 = n/6: 



"2(7) = ^cos 2 37. 



(2.14) 



Making the Taylor expansions around the minimum value 
of the gamma potential: 

9 9 81 / tt\2 „ / 7T\2 
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the equation for the variable 7 becomes 
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Using again in (|2.16[) the approximation 
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and making the change of variable y = 3j, we obtain 
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Eq. (|2.18|) is just the well known Mathieu equation. 
Using the expression for the characteristic value a, 
Ea. (|2.19[ ). of the Mathieu equation one can find the ex- 
pression for the excitation energy of the 7 equation 

fe2 1 
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The total energy for the system is obtained by adding 
the energies given by the equations (12.10[) and (|2.20l) : 

E{n p ,n 7 ,L,R) = E + Ep{np,L) + E 1 {n 1 ,L,R) (2.21) 

The excitation energies depend on four quantum num- 
bers, rip, n~f, L, R, and five parameters H 2 /2B, a, b, 

The quantum numbers defining the ground, beta and 
gamma bands are as follows: 
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(2.22) 

The wave function describing the whole system is: 
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where the factors -/V^ ^ and Ni,^ n denote the norms 
of the partial wave functions. 

Note that when the matrix elements with the wave 
functions depending on /3 and 7 respectively are calcu- 
lated, the integration over the /3 is performed with the 



measure f3 4 dft, while that over the 7 with the measure 
I sin 37^7. These measures are characterizing the (/?,7) 
space within the liquid drop model. These wave functions 
are further used to calculate the reduced E2 transition 
probabilities. 

In our approach the quadrupole transition operator is 
defined as: 
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The argument 7 — 27r/3 of the trigonometric functions 
is justified by the fact that it defines the axes 1 of the 
principal inertial ellipsoid. Indeed, the transformation 
from the laboratory to the intrinsic frame is a rotation 
defined by the matrix D^ R where the quantum numbers 
M and R are eigenvalues of the operator Q\, 

The reduced E2 transition probabilities are defined as: 



B(E2, J, -> Jf) 



l(^l|rf 2) ||J/)| 2 . 



(2.25) 



where the Rose's convention [24| was used for the reduced 
matrix elements. 

Summarizing, the SMA formalism uses a sextic oscil- 
lator potential with a centrifugal term for the j3 and a 
Mathieu equation for the 7 variable. These equations 
provide for the total energy given by Eq. (|2.2ip a com- 
pact form. The wave functions obtained by solving the 
quoted equations together with the transition operator of 
Eq. (|2.24p . are used to calculate the electric quadrupole 
transition probabilities. 

There are sever al g roups which studied the 7 soft nuclei 
around 7° = f [15T - [l8l |. The quoted approaches differ 
from the present formalism by the equations used for the 
description of /3 and 7 coordinates. 

Since the results of the SMA formalism will be com- 
pared with those obtained by CSM, in what follows we 
shall briefly present the main ingredients of the latter 
approach. 



III. COHERENT STATE MODEL (CSM) 

CSM defines 23] first a restricted collective space 
whose vectors are model states of ground, /3 and 7 bands. 
In choosing these states we were guided by some experi- 
mental information which results in formulating a set of 
criteria to be fulfilled by the searched states. 

All these restrictions required are fulfilled by the fol- 
lowing set of three deformed quadrupole boson states: 
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where the excitation operators for /? and 7 bands are 
defined by: 

nl = (6 w )o + J| ( 5t 6 t )o _^. (3 . 2) 



From the three deformed states one generates through 
projection, three sets of mutually orthogonal states 



where Pmk denotes the projection operator: 
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and Nj the normalization factors and D J MK the rotation 
matrix elements. It was proved that the deformed and 
projected states contain the salient features of the major 
collective bands. Since we attempt to set up a very simple 
model we relay on the experimental feature saying that 
the f3 band is largely decoupled from the ground as well 
as from the 7 bands and choose a model Hamiltonian 
whose matrix elements between beta states and states 
belonging either to the ground or to the gamma band 
are all equal to zero. The simplest Hamiltonian obeying 
this restriction is 



H 



A^N + 5nJ,fi/jO + A 2 J 2 + A 3 fijfi^, (3.5) 



where N is the boson number, J 2 -angular momentum 
squared and tip, denotes: 
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Higher order terms in boson operators can be added 
to the Hamiltonian H without altering the decoupling 
condition for the beta band. An example of this kind is 
the correction: 



AH = A^tiptij, + h.c.) + A 5 tif,tij,. 
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The energies for beta band as well as for the gamma 
band states of odd angular momentum are described as 
average values of H (|3.5[) . or H + AH on (flj M and <pj M 
(J-odd), respectively. As for the energies for the ground 
band and those of gamma band states with even angu- 
lar momentum, they are obtained by diagonalizing a 2x2 
matrix for each J. 

The quadrupole transition operator is considered to 
be a sum of a linear term in bosons and one which is 
quadratic in the quadrupole bosons: 

Q2, = 51 (4p + i-TK-,) + qiiftb^ + {bb) 2 ,). (3.8) 



The form of the anharmonic component of Q 2 ^ is justified 
by the fact that this is the lowest order boson term which 
may connect the states from beta and ground bands in 
the vibrational limit, i.e. d-small. 

Using the Rose convention [24j , the reduced probability 
for the E2 transition J" 1 



can be expressed as: 



B(E2;J^Jp=[(J+\\Q 2 \\Jp 
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Three specific features of CSM are worth to be men- 
tioned: 

a) The model states are generated through projection 
from a coherent state and two excitations of that through 
simple polynomial boson operators. Thus, it is expected 
that the projected states may account for the semiclas- 
sical behavior of the nuclear system staying in a state of 
high spin. 

b) The states are infinite series of bosons and thus 
highly deformed states can be described. 

c) The model Hamiltonian is not commuting with the 
boson number operator and because of this property a 
basis generated from a coherent state is expected to be 
most suitable. 

The CSM has been successfully applied to several nu- 
clei exhibiting various equilibrium shapes which accord- 
ing to the IBA (Interacting Boson Approximation) classi- 
fication, exhibit the SO(6), SU(5) and SU(3) symmetries, 
respectively. Several improvements of CSM has been pro- 
posed by considering additional degrees of freedom like 
isospin [25] . quasiparticle [26j j or collective octupole co- 
ordinates^. [28j . CSM has been also used to describe 
some nonaxial nuclei [29j and the results were compared 
with those obtained with the Rotation- Vibration Model 

3C I . A review of the CSM achievements is found in Ref. 

31|. 



IV. SMA FORMALISM OBTAINED BY 
QUANTIZING THE CLASSICAL CSM 
EQUATIONS. 

In our previous publication on this subject [l9| we were 
noting that the two approaches SMA and CSM describe 
the data of the considered nuclei, equally well. This 
amazing feature raised the question why that happens? 
Actually, here we aim at answering the mentioned ques- 
tion. In brief we shall prove that, indeed, SMA can be an- 
alytically derived from the CSM formalism. The source of 
our inspiration was the result from Ref. [32| showing that 
a generalized Davidson potential [33[ can be obtained by 
a semiclassical treatment of a fourth order boson Hamil- 
tonian. 

For this purpose we shall study the classical properties 
emerging from CSM, by dequantizing the specific boson 
Hamiltonian used by CSM, considering its average with 
the coherent state: 



\ip) = exp z bl + z 2 b\ + Z- 2 b'_ 
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where Zk,z£ with k = 0, ±2 are complex numbers de- 
pending on time. As usual, |0) denotes the vacuum state 
for the quadrupole bosons. The basic property of this 
function is comprised by the equation: 



(4.2) 



The classical Hamilton function associated to the CSM's 
model Hamiltonian is: 



U = (i/>m) = 2(11A 1 +3A 2 ) (|z | 2 + M 2 + k- 2 | 2 ) 
+ A x (2z 2 z*_ 2 + zf - d 2 ) (2z 2 z_ 2 + zl - d 2 ) + 

+ ^ [2 (6z*z* 2 zt 2 - z* 3 ) + M (2z*z*_ 2 + z* 2 ) - d 3 } 

x [2 (6z z 2 z_ 2 - zl) + 3d (2z 2 z_ 2 + z 2 ) - d 3 } . (4.3) 

The equations of motion described by the classical coor- 
dinates zq, z±2 and their complex conjugates Zq,z± 2 are 
obtained from the variational principle of the minimum 
action: 



5 j\^\H-ih^p\i-),ll' = n. 
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The action variation can be written in a compact form: 
6(ip\iH-^-\ip) = z 5zq-z Sz q 

. * 

+ z 2 Sz 2 ~ z 2 Sz 2 

+ z_ 2 5z*_ 2 — z_ 2 8z^ 2 . (4.5) 

where the symbol " ■ " stands for the time derivative. 
Finally, the result for the classical equations is: 
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These equations support the interpretation of z^, with 
k = 0, ±2, as classical phase space coordinates and of z%. 
as the corresponding conjugate momenta. With the com- 
plex coordinates we may define the canonical conjugate 
coordinates: 
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V2 V2 
which, evidently, obey the equations: 

{Qn,Po} = 1, {Q±,P±} = 1, 

{Qk,H} - Q k , {P k ,H}=P k . (4.8) 

For two given functions /, g defined on the phase space, 
their Poisson bracket is denoted by {/, g} and defined as: 



{/,<?} = 



v ( df dg df dg 
^ \dQ k 8P k dP k 8Q k 
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In terms of the canonical coordinates Q and P the clas- 
sical energy function is given in Appendix A. 

In what follows we shall study the Hamilton function 
H in the subspace defined by z 2 — z- 2 , where we use the 
canonical coordinates defined by: 



Qo = Qo, Po = Po, 
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Q 2 + Q- 
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These coordinates are related with the real and imagi- 
nary part of the complex variable Zk, by the following 
equations: 

q = V2u , po = V2v Q , q 2 = 2u 2 , p 2 = 2v 2 , (4.11) 
uq = Re zq, vq = Im zq, u 2 — Re z 2 , v 2 — Im z 2 . 

The restriction to the mentioned subspace is justified 
by the following considerations. It is well known the 
fact that due to the overcomplete property the coherent 
states may be used to construct basis functions describ- 
ing specific irreducible representations. In particular the 
coherent state used in this paper can be used to project 
out the basis \NvaJM) where the quantum numbers are 
the boson number N, the seniority v, the missing quan- 
tum number a, the angular momentum J and its projec- 
tion on z-axis M (34J. In Ref.[34j| it was shown that a 
two parameters coherent state is sufficient to be able to 
project out the entire basis mentioned above. Therefore 
the properties of the GSM Hamiltonian in the whole bo- 
son space can be described in the reduced phase subspace 
mentioned above. 

In the expression of H we adopt the approximation 
and coordinate transformations: i) we neglect the terms 
non-quadratic in momenta as well as the therms coupling 
the coordinate with momenta; ii) we take care of the 
restriction z 2 — z-2 by introducing the new canonical 
conjugate coordinates (qo;po) and (g2;P2); hi) for the new 
coordinates we use the polar coordinates: 



qo = r cos 7, q 2 = r sm 7 



(4.12) 



In this way H. is a sum of the kinetic and potential energy 
terms: 

H = (lUi + 3A 2 + A,d 2 + ^d 4 A 3 )( P 2 + p\) + V(r, 7), 
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In order to obtain a separable equation for the vari- 
ables r and 7 we approximate V(r, 7) by a sum of two 
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potentials one depending only on r, Vi(r), and the other 
one only on 7, ^2(7). In the terms of Vi(r), the factors 
depending on 7 are considered in the minimum point of 
^2(7) (which is 7r/6), while in ^2(7) the factors depend- 
ing on r are considered in the minimum point of V\ (r), 
denoted by tq. The approximated potential will be de- 
noted by U(r 7 7). 
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In Figs. 1,2, panels (a), (c), (e) we present the con- 
tour plot and two sections of the approximate potential 
U(r, 7). The two sections are obtained by fixing one of 
the two variables in the potential minimum point. In 
the panels (b), (d) and (e) similar plots for the exact 
potential V(r, 7) are given. Figures 1 and 2 correspond 
to different sets of the CSM parameters, i.e. d,A x ,A2 
and A 3 . Comparing the curves from left and right panels 
of each figure we conclude that the effect of approxima- 
tions yielding the separated form of the potentials is quite 
small. The only visible effect is on the orientation of the 
symmetry axis of the equipotential curves. Indeed, the 
contour plots for exact and approximated potentials have 
different symmetry axes although they surround the same 
minimum point. Changing slightly the CSM parameters 
we note that indeed the curves from Fig. 2 are quite 
close to the corresponding curves from Fig.l. Thus, we 
may say that indeed, the curves are stable against small 
changes of the set of parameters determining the classical 
r and 7 potential. 

The classical Hamilton function becomes: 



u 



(HAi + 3A 2 + Aid 2 + —d 4 A 3 ){p 2 +p 2 2 ) + U(r, 7). 



(4.15) 

This can be quantized by replacing the sum of the mo- 
menta squared by the Laplace operator written in polar 
coordinates: 
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Vi(r) + V 2 ('y). 



It is convenient to introduce the notation: 
T = (lLAi + 3A 2 + A x d 2 + ^-d A A 3 ) 
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The Schrddinger equation: 



for the trial function: 

*(r, 7 ) = Vi(r)e^2(7), 
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FIG. 1: (Color online) The contour plot for the approximated 
potential U(r, 7), given in (4.14), is presented in the variables 
r and 7 (a). The parameters involved are d = 3.5, Ai = 
W.65keV, A 2 = 15.U6keV, and A 3 = 150/fceV. The min- 
imum is reached in (ro,7o) = (2.6, 7r/6). The equipotential 
curves are separated by an amount of 80 keV. The mini- 
mum value of the potential is U m in = — 3305.25fceV. Two 
sections 7 = f and ro = 2.6 of the potential U(r, 7) are 
presented in panels (c) and (e), respectively. Similar plots 
but for the exact potential from (4.13) are given in the pan- 
els (b), (d) and (f), respectively. Parameters were kept the 
same as for the approximated potential. The minimum value 
of the potential is V m m = — 3280.31fceV and is reached at 
(>"o,7o) = (2.7, 7r/6)..The equipotential curves are separated 
by an amount of 50 keV. 



is separated 
'I d 



T ™ 



2 



r dr dr 2 r 
UK d ^ 1 d 2 
r 2 dj 2 



dj 



+ V 1 (r) %j) X (r) = E^ipi{r), 

+ ^2(7) Ml) =E^M7)- 
(4.20) 

In what follows we shall show that the first equation 
(|4.20p leads to a Schrddinger equation for a sextic po- 
tential plus a centrifugal term, while the second equation 
provides a differential equation obeyed by the Mathieu 
function. 

Dividing both sides of the first Eq. (|4.20l) by T and 
denoting by 

u 1 {r)=V 1 {r)/F, e r = E^/T (4.21) 
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FIG. 2: (Color online) The contour plot for the approximated 
potential U(r,"f) given by Eq. (4.14) is represented as func- 
tion of the variables r and 7 (a). The parameters involved are 
d = 3.22, Ai = 12.95fceF, A 2 = 7.9fceV, and A 3 = 240fceV. 
The minimum is reached in (ro,7o) = (2.4, £). The equipo- 
tential curves surrounding the minimum value are separated 
by an amount of 70 keV. The minimum value of the po- 
tential is Umin = — 3276.43/ceV. Two sections 7 = ~ and 
ro = 2.4 of the potential U(r, 7) are presented in the panels 
(c) and (e), respectively. Similar plots but for the exact po- 
tential from (4.13) are given in the panels (b), (d) and (f), 
respectively. Parameters were kept the same as for the ap- 
proximated potential. The minimum value of the potential 
is Vmin = — 3250.92fceV r and is reached at (ro,7o) = (2.5, -|). 
The equipotential curves are separated by an amount of 50 
keV. 



the equation in the variable r becomes: 

Mr)=£rMr)- ( 4 - 22 ) 



d 2 1 d K 2 



dr 2 r dr r 
Changing the function 

i>i (r) = r~i4>{r), 
the equation for the new function is 



d 2 

Or 



K 2 -\ 



ux(r) 



(4.23) 



(4.24) 



which is nothing else but the Schrodinger equation for a 
sextic potential plus a centrifugal term. 



Now let us turn our attention to the second equation 
from (|4.20[) . Multiplying it with r^/T and denoting by: 



_ ^3 8 
M ~ 280 J" r °' £l 



]E<?>_ 
T 



(4.25) 



d 

2iK— — h fx cos 67 

07 



one obtains: 
" 

<9 7 2 

With the change of function: 

Ml) = e"^M(3 7 ) 

we obtain: 

d 2 

dl 2 



Ml) = ^Ml)- (4-26) 



(4.27) 



7—2 + e-y + K — fj, cos 67 



M(3 7 ) = 0. 



(4.28) 



Changing now the variable 7 to y — 37, the equation for 
the Mathieu function is readily obtained: 



dy 



+ a — 2q cos 2y 



M(y) = 0, 



where the following notations have been used: 



(4.29) 



(4.30) 



Before closing this section we would like to comment 
on the relation of the variable r and the dynamic nuclear 
deformation B. 

Aiming at this goal let us consider the canonical trans- 
formation relating the quadrupole conjugate coordinates 
and the boson operators: 



«2 M 



7T2, 



75 (bl + (-Vh,^) , 



kV2 

- — (V (-Y 



r )2,i 



(4.31) 



Note that the canonical transformation from above is de- 
termined up to a multiplicative factor k. Averaging these 
equations with the coherent state ip (5.1), one obtains 
that the coordinates Q M and P M introduced above are 
related with the quadrupole operators by: 



(ifj\ka 2 ^\ip), 



(4.32) 



Identifying the averages of with the coherent state, 
with the quadrupole coordinates in the intrinsic reference 
frame we obtain (35| : 

kB 

Qo = ka 2 Q = kBcos-f, Q±2 = — psin7. (4.33) 

v2 

In the restricted classical phase space the canonical co- 
ordinates are: 



go = fc/3cos7, »2 = kB sin 7. 



(4.34) 
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From here it results: 



r = k/3 



(4.35) 



Using this simple relation in connection with the differ- 
ential equation in r, one obtains the Shrodinger equation 
for sextic potential plus centrifugal term in the variable 



'd/3 2 



K 2 



2 



k 2 ui{k/3) 



= Ef)<f>(kp) 



(4.36) 

Concluding this section, we may say that while in the 
previous paper [l9| the sextic potential for f3 and the 
7 potential yielding the equation for the Mathieu func- 
tion were introduced by an ad-hoc choice, here they are 
derived in a natural manner from the CSM formalism. 
Moreover the variable separation is based on two ap- 
proximations suggested by the classical picture: i) the 
non-quadratic terms in momenta are ignored and ii) the 
coupling of coordinates and momenta are vanishing due 
to the local character of the classical phenomenological 
forces. 



V. NUMERICAL RESULTS 

The formalisms SMA and CSM presented in the pre- 
vious sections have been applied for calculating the ex- 
citation energies and the available B(E2) values for two 
isotopes: 180 Hf, 182 W. We start with the excitation en- 
ergy analysis. As shown above the total energy provided 
by SMA depends on five parameters: h 2 /2B, a, b, jpzji A*- 
These have been fixed by fitting the excitation energies 
using the least square procedure. The results are given 
in Table I. Concerning CSM, the parameters determin- 
ing the energies are: d, A\, A2, A3. They were fixed as 
follows. We cycled d within a large interval with a small 
step. For each d we determined A\ and A2 by fitting 
the energies of two states, one belonging to the ground 
band and one from the gamma band. A3 was obtained 
by fitting one level energy from the beta band. Then we 
choose that d which yields an overall good fit. The fitted 
parameters are given in Table I. 

The results obtained with the two approaches, SMA 
and CSM, are compared, in Fig. 3 and Fig. 4, with the 
corresponding experimental data. 

By inspecting the experimental data we may see 
whether the signature for triaxiality shows up. Indeed, 
it is well known the fact that the most distinctive signa- 
ture of the triaxial rigid rotor is the equation relating the 
energies of three particular states [36p[38| : 



E 2 + + E 2 + = E 3 + . 



(5.1) 



Actually this equation is only approximately obeyed. De- 
noting by AE the modulus of the difference between the 
left and right hand side of the mentioned relation, the 
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g band 

+ 2248 .2272. .2273. 
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f3 band 
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1100 



6* 
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1084 1083 
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y band 

1568 1557 1544 

1+ 1414 1409 1401 
, + _1292_ 1291 1286 

t [183 1200 W + U6Q 1183 1188 
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180 
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2 + . 
+ . 



99 93 




93 



SMA Exp. CSM 



FIG. 3: Excitation energies for ground, beta and gamma 
bands in 180 Hf, obtained with SMA and CSM formalism re- 
spectively, are compared with the corresponding experimental 
data taken from Ref. 391. 



experimental data lead to the values: 

AE = 2keV; lOkeV. (5.2) 

for 180 Hf and 182 W, respectively. Clearly, these devia- 
tions suggest that the nuclei considered in the present 
paper are close to an ideal triaxial rotor. As a matter of 





180 Hf 


182 W 


£[keV] 


0.401 


0.476888 


a 


16212.86 


11508.56 


b 


-44 


95 


W) 


3.14 


3.457 


ft 


20892. 


12772.6 


h 


177.93[W.u] 1/2 


158.48 [W.u] 1/2 


h 


4630.24[W.u.] 1/2 


4736.037 [W.u.] 1/2 


d 


3.5 


3.22 


Ai[keV] 


21.17 


21.54 


A 2 [keV] 


8.15 


7.47 


A 3 [keV] 


-10.85 


-12.29 


Qi 


3.739[W.u.] 1/2 


3.756[W.u.] 1/2 


92 


-0.125IW.U.] 1 / 2 


-0.175W.U.] 1 / 2 



TABLE I: The parameters h 2 /2B, a, b, 



7TTT 



fi involved in 



the energy expression provided by SMA (|2.21 [1 . are given for 
180 Hf and 182 W. Also we give the values for the parameters ti 
and t 2 defining the transition operator used by SMA H2.24|) . 
On the last six rows we give the parameters determining the 
CSM excitation energies, d, Ax, A<z, A3, and the specific E2 
transition operator i.e., q\ and 92- 
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FIG. 4: Excitation energies for ground, beta and gamma 
bands in 182 W, obtained with SMA and CSM formalism re- 
spectively, are compared with the corresponding experimental 
data taken from Ref . [4011 . 



fact this is the experimental feature which inspired us to 
take the 7 = 30° as reference picture. 

Let us discuss now the results concerning the transi- 
tion probabilities. The SMA made use of a anharmonic 
transition operator written in the intrinsic frame of ref- 
erence (|2.24[) . while CSM employs a second order boson 
operator in the laboratory frame (|3.8I) . In both cases the 
operators involve two parameters: t\ and t 2 for SMA and 
qi and q 2 for CSM. These parameters have been fixed by 
fitting two particular transitions for each nucleus. The 
fitted parameters are given in Table I. Results of our cal- 
culations and the corresponding experimental data are 
given in Table II for 180 Hf and Table III for 182 W. 

We remark that both formalisms describe reasonable 
well the reduced transition probabilities. It is worth not- 
ing the high accuracy of the CSM description of both 
excitation energies and transition probabilities for 180 Hf. 

Before closing this section we would like to spend few 



TABLE II: B(E2) values for some ground to ground and 
gamma to ground E2 transitions in 180 Hf. J+ and J+ denote 
the angular momenta of the initial and final states, respec- 
tively. Experimental data were taken from Ref. [391] 



words about the other descriptions of the triaxiality fea- 
tures. We have not to forget however the fact that the 
pioneering paper for phase transition from gamma stable 
to gamma unstable nuclei with an analytical description 
of the critical point, and the departure from axial sym- 
metry is that of Jean- Wilets [4l| . However an intensive 
study of the subject started in the beginning of the last 
decade, since the context of symmetries was much de- 
veloped and on the other hand relevant data have been 
accumulated. Triaxiality has been investigated within 
the IBA formalism being related to various effects. Thus 
including higher order terms the triaxiality of 190 > 192 Os 
has been studied in Ref. [42j]. Including the g-boson, 
in a recent study (43| no shape/phase transition towards 
stable triaxial shapes has been found. The phase dia- 
gram of IBA-2 (which distinguishes protons from neu- 
trons) including triaxial shapes, has been constructed in 
Refs. gilig. 

Several authors treated the gamma soft nuclei around 
7° = § [T5MT8I ] . However, their equations for beta as well 
for gamma variables are different from those proposed in 
the present paper. The beta potential is either an infi- 
nite square well [15J or a Coulomb or a Kratzer potential 
[16l - tl8l ]. Recently [43], the Davidson potential was used 
in relation to triaxial nuclei. Concerning 7 all quoted de- 
scriptions use an oscillator potential. The sextic poten- 
tial was previously used in Refs. [48], [49| but only for few 
low lying states. Again the description of 7 is different. 
Triaxiality has recently been studied in the framework of 
the algebraic collective model [50[ , and the onset of rigid 
triaxial deformation has been considered (5l| . 

The structure of the projected states in the intrinsic 
variable is relevant for the subject under consideration. 
In Ref. [52j we studied the probability distributions of 
the 7 and the /3 variables corresponding to various states 
from the ground, beta and gamma bands. For the state 
0~p the 7 probability has a minimum in ir/6 and two max- 
ima in 7 = and 7 = 7t/3, respectively. Increasing the 
spin in the (3 band, the picture is changed. For example 
1+ 



for 10+ the probability distribution of 7 has a minimum 
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► 16+ 
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315 




► 0+ 
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200 
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>o+ 


1.3 


0.9 


0.008 






10.5 


1.7 


0.021 




>o+ 


0.2 


3.4 


3.4 




►2+ 


8.5 


6.74 


6.27 




►4+ 


0.0 


0.034 


0.51 




>2 + 


0.1 


2.4 


1.36 




►4+ 


1.7 


10.4 


7.60 



TABLE III: B(E2) values for some ground to ground, beta to 
ground and gamma to ground E2 transitions in 182 W. J+ and 
J+ denote the angular momenta of the initial and final states, 
respectively. Experimental data were taken from Ref. [||| • 



in 7r/6 and two maxima, one for 7 = 7r/12 and one for 
7 = 7r/4. The 7 band exhibits an oposite feature, namely 
the head state has a maximum in 7r/6, while in a high 
state, like J 71 " = 10+, two preeminent maxima for 7 = 
and 7 = 7r/3 are observed. In the ground state 0+ one 
meet a situation which is specific to gamma unstable nu- 
clei, i.e. the distribution is almost constant. Increasing 
the spin one obtains for J T = 10+ a maximum at 35°. 
In beta variable the state 0^" has a bimodal structure, 
the head state of gamma band has a deformed maxi- 
mum while the ground state a spherical maximum. If 
the factor /3 , entering the measure of integration over 
the dynamical deformation /3, then the distributions are 
peaked on various deformations, depending on the band 
to which the state belongs and of course on the angular 
momentum. The resuts mentioned above are consistent 
with the predictions of the semiphenomenological model 
proposed by Kumar and Baranger [53| . 

Recall the fact that the ground band states with J > 2 
are mixed with the states of the same angular momentum 
from gamma band. Therefore one expects to have gamma 
distributions with a maximum value at 7 = 7r/6, at least 
for low lying states from the ground and gamma bands, 
while in the beta band two maxima at about 7r/12 and 
7r/4 may show up. 

The transition from a near spherical shape to 
a triaxial shape is accompanied by the change 
of the staggering from 2+, (3+, 4+), (5+ 6+), ... to 
(2+,3+),(4+,5+),(6+,7+),..., In Ref. || we showed 
that whithin CSM the first clustering is remnant of the 



doublet degeneracies of the vibrational limit of the model, 
while the second staggering which is typical for rigid tri- 
axial rotor, is obtained in the asymmptotic region of de- 
formation. Also, such a transition is reflected in the be- 
haviour of the E2 transitions in the gamma band. 

We remember that SMA was initially [19( obtained 
from the Bohr-Mottelson Hamiltonian by separating the 
variables j3 and 7 and supplementing the result with 
new potentials for beta and gamma respectively. Here 
SMA is derived from CSM based on a similar procedure 
of the variable separation. Of course it would be nice 
if the final quasiexactly separable approach keeps track 
of the formalism from which it emerges. This problem 
was studied in great details by Caprio [54] in connec- 
tion with the X(5) Hamiltonian comparing the results of 
approximate separated Hamiltonian with the exact ones 
obtained by diagonalizing the initial Hamiltonian in a 
large basis constructed from a five dimensional spherical 
harmonics [55T4591 ] . The conclusion was that replacing 
1//3 2 by the "rigid" value 1//3q is a valid approximation 
for nearly 7 soft potential, while the "small" angle ap- 
proximation for 7 is good for large 7 stiffnesses. Both 
approximations work in the overlapping region of the two 
7 intervals. 

Within the SMA the situation is completely differ- 
ent. The only term approximated in the 7 equation is 
1 / sin 2 37 originating from the 7 kinetic energy term. 
However, this is not expanded around 7 = which is 
related to the 7- small picture but around 7r/6 where 
the mentioned term is minimum. In Ref. [19| we studied 
the potentials in 7 for several nuclei and noticed that 
they are characterized by a large stiffness. Concerning 
the approximation in /3, in the quoted paper we plotted 
the functions of different states from the three bands. 
For all considered cases, the functions are highly local- 
ized on /3o- There is an esential difference between our 
case and that of Ref. [H^. Indeed, there the Hamiltonian 
which is separated is harmonic and the approximations 
concern the part of centrifugal term 1//3 2 which is dis- 
tributed to the factor expressing the 7 kinetic energy and 
the quadratic Taylor expansion around 7 = giving the 
fact that this is a singularity point for the 7 potential. 
Note that in our descrition the Hamiltonian is highly an- 
harmonic and moreover the 7 expansion is made around 
a stability point for the 7 potential. Moreover the /3 
potential has a deep deformed minimum which favor a 
small contribution of fluctuations which might shift the 
average (/3 2 ) from (3%. 

Concerning the SMA derivation from the CSM, we 
don't have a centrifugal term in r which is coupled to 
the terms from the 7 kinetic energy, and therefore the 
trouble generated by replacing r with r$ in a centrifu- 
gal like term is absent. However, the separation is ulti- 
mately obtained by replacing 7 by the minimum point 
coordinate 70 in some terms, and r with the minimum 
coordinate ro in other terms. Such an aproximation is 
justified for situations when the eigenfunctions of the 7 
and r Hamiltonian respectively, are well centered on the 
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7o and r$. This picture might be reached indeed, since 
as shown in Figs. 1, 2, the potentials depths are large. 
Actually, a confident measure of the separation approx- 
imation quality is the comparison of the approximated 
and non-approximated potentials shown in Figs 1, 2 for 
two sets of parameters. As seen there, the potentials do 
not differ from each other too much as to induce a sig- 
nificant difference of the corresponding wave function. 



12 



VI. CONCLUSIONS 

Here we summarize the main results obtained in the 
previous sections. The formalism, conventionally called 
Sextic and Mathieu Approach (SMA), proposed in a pre- 
vious publication for the description of the triaxial nuclei 
is applied to describe the spectra and transition proba- 
bilities for another two nuclei, 180 Hf and 182 W. Results of 
our calculations are in good agreement with the available 
data. 

The main result of the present study is a natural 
derivation of the SMA as a limiting case of the CSM. The 
model Hamiltonian of CSM together with a three param- 
eters quadrupole coherent state is used in a time depen- 
dent variational principle to derive the classical equations 
of motion in the classical phase space coordinates. In the 
classical energy function one neglects the terms which are 
non-quadratic in momenta as well as those which couple 
coordinates and momenta. Reducing the space to the 
subspace generated by two pairs of conjugate coordinates 
and then quantizing the classical Hamilton function one 
arrives at a separable form for the associated Schrodinger 
equation, one being the equation for a sextic oscillator in 
ft and the other one a differential equation for 7 obeyed 
by the Mathieu function. The reduced space corresponds 
to a two parameters coherent state from which the entire 
boson space basis \NvaJM) can be projected out. This 
justifies the restriction since in this way one accounts for 
the classical properties of the CSM model Hamiltonian 
in the whole boson space. 

The hnal conclusions are: 

• The sextic potential for the /3 variable and the po- 
tential in 7, introduced in our previous paper on 
this subject [19j, based on pragmatic grounds, gets 
a theoretical support. 

• The success of CSM in explaining the data in a re- 
alistic way lets us conclude that, indeed, this model 
is able to describe the triaxial nuclei. 

• According to the results of the present paper, the 
SMA represents the image of CSM in the intrinsic 
reference frame. 

• We stress, again, that in describing the near vi- 
brational, well deformed with axial symmetry, tri- 
axial or gamma unstable nuclei one uses a sole 
model Hamiltonian and a sole set of model func- 
tions for the states belonging to the ground, beta 
and gamma bands. 

• This makes us hope that the phase transition from 



prolate to oblate shape through the triaxial shape 
can be described by CSM. 

This feature will be however presented, soon, in another 
publication. 
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H = (lUi + 3v4 2 )(P 2 + 2P_ 2 P 2 + Q 2 + 2Q_ 2 Q 2 ) 

fj) 

4 



Ai(d 4 + d 2 P 2 + -f + 2d 2 P_ 2 P 2 + P- 2 P 2 P 2 



P 2 2 P| + P\Q\ + 2P_ 2 P Q_ 2 Q - d 2 Ql 
ip 2 Q 2 - P-2P2QI + % - 2d 2 Q- 2 Q 2 - P0Q-2Q2 
2P0P2Q0Q2 + Q-2QIQ2 + PlQ\ + QUQl) 
^ + 3d 4 P 2 + \d 2 P± + % + 6d 4 P_ 2 P 2 
9d 2 P_ 2 P 2 P 2 - 6P_ 2 P 4 P 2 + M 2 P 2 2 P| 
18P 2 2 P 2 P 2 2 + 6V2d 3 P_ 2 P Q_ 2 + 12V2dP_ 2 PoQ-2 
9d 2 P 2 2 Q 2 _ 2 + 18P 2 2 P%Q 2 _ 2 - 3V2d 3 P£Qo 

+ 6V2d 3 P_ 2 P 2 Q - 18V2dP_ 2 P 2 P 2 Q 

V 2 

18V2dP 2 2 P 2 Q + 18rf 2 P_ 2 P Q_ 2 Qo - 12P_ 2 P 3 Q_ 2 Q 
18V2dP 2 2 Q 2 _ 2 Q - id i Ql + ^d 2 P 2 Ql 



^Q 2 - 9d 2 P^ 2 P 2 Q 2 + 18P 2 2 P 2 Ql 
\8P 2 2 Q 2 _ 2 Ql + V2d 3 Q 3 - 3V2dP 2 Q 3 - 6\/2dP_ 2 P 2 Q 3 
12P_ 2 P Q_ 2 Q 3 + ^d 2 Q 4 a + ^P 2 Q 4 + 6P_ 2 P 2 Q 4 

^3 + ^ + 6v / 2d 3 P P 2 Q 2 + 12V2dP 3 P 2 Q 2 
V 2 2 

6d 4 Q_ 2 Q 2 - M 2 P 2 Q_ 2 Q 2 + 6P 4 Q_ 2 Q 2 
18d 2 P P 2 QoQ 2 - 12P 3 P 2 QoQ 2 - 6V2d 3 Q_ 2 QoQ 2 
18V2dP 2 Q_ 2 Q Q 2 + 9d 2 Q_ 2 Q 2 Q 2 - 12P P 2 Q 3 Q 2 
6%/2dQ_ 2 Q 3 Q 2 - 6Q_ 2 g^Q 2 + 9d 2 P 2 Q 2 
18P 2 P 2 Q 2 2 + 9d 2 Q 2 _ 2 Q 2 2 + 18P 2 Q 2 _ 2 Q 2 2 + l8V2dP 2 Q Q 2 2 
18V2dQ 2 _ 2 QoQ 2 2 + \8PlQlQ\ + 18Q 2 _ 2 Q 2 Q 2 2 ) . (A.l) 
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